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I.  Introduction 

In  the  Soviet  Union  the  method  of  artificially  freezing  ground 
is  being  used  more  and  more  every  year  in  the  construction  industry. 

Under  complex  hydrogeological  conditions  this  method  makes  it  possible 
to  create  water-impermeable  retaining  walls,  cofferdams  and  arches, 
and  ground  can  be  reinforced.  The  area  in  which  ground  freezing 
can  be  applied  includes  all  the  basic  branches  of  the  construction 
industry.  Some  idea  of  the  work  which  has  been  carried  out  can  be 
gained  from  the  fact  that  in  the  last  20  years  in  our  country  more 
than  l.S  million  cubic  meters  of  ground  have  been  frozen  [1]. 

It  is  necessary,  however,  to  admit  the  fact  that,  despite  the 
tremendous  volume  of  the  work  which  has  been  carried  out,  the  level 
of  our  knowledge  of  freezing  technology  has  not  yet  reached  the 
necessary  level.  This  fact  became  especially  evident  when  the 
freezing  method  was  used  on  several  hydrotechnical  construction 
projects  where  water-permeable  soils  were  subjected  to  freezing. 

Recent  studies  carried  out  by  the  Ice  Thermics  Laboratory  of 
the  B.  Ye.  Vedeneyev  All-Union  Scientific  Research  Institute  of 
Hydraulic  Engineering  have  shown  that  the  methods  used  by  design 
organizations  to  make  thermal  calculations  of  ground  freezing  and 
primarily  of  the  freezing  of  water-permeable  soils  contained 
significant  flaws.  In  this  process  it  was  found  that  one  of  the 
basic  shortcomings  is  the  extremely  approximate  allowance  which  is 
made  for  the  thermal  influx  from  the  unfrozen  ground  to  the  surface 
of  the  frozen  mass.  When  we  examine  the  problem  of  freezing  ground 
in  which  there  is  no  water  motion,  this  drawback  means  that  the  actual 
freezing  time  will  not  coincide  with  the  projected  data;  however, 
freezing  under  these  conditions  is  usually  successful  due  to  the 
fact  that  the  intensity  of  the  heat  flow  from  the  outside  decreases 
with  time.  Therefore  the  existing  methods  of  calculating  the  freezing 
rate  of  nonwater-permeable  soils  are  adequate  for  engineering 
practice  to  a certain  extent. 

However,  the  picture  is  different  when  we  discuss  the  freezing  of 
water-perneable  ground.  In  this  case  the  thermal  flux  from  the  outside 
does  not  decrease  with  time,  and  as  a result  of  this  it  may  be 
found  that  a design  solution  based  on  calculation  results  in  which 
the  thermal  flux  from  the  outside  is  considerably  lower  than  it  is 
in  reality  will  be  proven  wrong  in  practice  (for  instance,  it  will  be 
impossible  to  make  the  frozen  ground  cylinders  join  or  the 
assigned  thickness  of  the  frozen  ground  walls  will  not  be  achieved). 
Therefore  it  is  natural  to  direct  special  attention  at  improving  the 
methods  of  making  thermal  calculations  of  the  freezing  of  water-permeable 
soils. 

In  1951  the  pages  of  the  "B.  Ye.  Vedeneyev  All-Union  Scientific 
Research  Institute  of  Hydraulic  Engineering"  published  a work  by 
B.  V.  Proskuryakov  which  provided  a new  method  for  thermally  calculating 
the  freezing  of  the  water-permeable  ground  by  a single  column  [2]. 
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A subsequent  experimental  check  made  by  the  Institute  of  this 
calculation  method  proved  that  the  author's  equations  were  correct 
and  simultaneously  demonstrated  the  great  imprecision  of  the  formulas 
which  had  been  used  up  to  that  time. 

At  the  present  time  it  is  generally  accepted  that  the  speed 
with  which  ground  is  frozen  by  a group  of  columns  before  the  frozen 
ground  cylinders  join  is  carried  out  without  allowance  for  the 
mutual  influence  (interference)  of  the  frozen  ground  cylinders,  i.e., 
the  freezing  of  the  ground  by  one  group  column  is  generally  equated 
with  the  freezing  of  ground  by  a single  column.  It  must  be  noted 
that  this  is  completely  wrong.  Thus,  for  instance,  in  reality  when 
freezing  nonwater-permeable  ground  a single  column  will  have  less 
favorable  conditions  than  a group  column  in  view  of  the  fact  that 
the  frozen  ground  cylinder  of  a group  column  is  partially  screened 
with  regard  to  heat  by  its  neighboring  cylinders  and  therefore  in  this 
case  the  thermal  influx  from  the  unfrozen  ground  will  be  less  than 
in  the  case  of  a single  frozen  ground  cylinder.  In  water-permeable 
ground  the  presence  of  neighboring  cylinders  also  exerts  an  influence. 

We  consider  it  very  desirable  that  in  the  near  future  the  problem  of 
the  mutual  influence  of  freezing  columns  both  in  water-permeable 
and  in  nonwater-permeable  soils  be  studied  in  detail.  In  the  meantime, 
however,  we  will  assume  that  the  calculations  of  ground  freezing 
rates  before  the  solidification  of  the  frozen  ground  cylinders 
will  be  carried  without  allowance  for  the  influence  of  the  above- 
-mentioned  column  interference  factor. 

In  this  article  we  shall  examine  the  problem  of  ground  freezing 
after  the  frozen  ground  cylinders  join.  In  doing  so  we  will  state  and 
solve  only  a single  partial  problem:  when  the  freezing  columns  are 
located  in  a single  straight  row  and  the  filtering  flow  after  the 
frozen  ground  cylinders  join  is  directed  along  the  frozen  ground 
wall . 

2.  Statement  of  the  Problem 

The  above-mentioned  problem  of  creating  a frozen  ground  wall 
can  be  formulated  as  follows. 

A row  of  freezing  columns  is  implanted  in  an  infinite,  homogeneous 
layer  of  water-permeable  ground  perpendicular  to  the  layer  (Figure  1). 
The  top  and  bottom  of  the  layer  are  two  mutually  parallel  surfaces  which 
are  impermeable  to  water  and  heat.  The  length  of  each  column  is 
equal  to  the  thickness  of  the  layer  h.  Columns  with  radius  rQ  are 

arranged  with  interval  F between  them.  The  temperature  of  the 
filtration  flow  at  a considerable  distance  from  the  column  is  equal 
to  T.  After  the  frozen  ground  cylinders  join,  the  filtration 
flow  is  directed  along  the  frozen  ground  wall.  The  speed  of  the 
filtration  flow  is  equal  to  v.  The  ground  freezes  at  temperature  t^. 

The  physical  characteristics  of  the  frozen  and  unfrozen  ground  are 
known.  It  is  necessary  to  find  a dependency  which  determines  the 
freezing  speed  of  the  ground  after  the  frozen  ground  cylinders 
join. 


Horizontal  cross-section  of  frozen 
layer  of  ground 


e O' 


Vertical  cross-section  of  frozen 
layer  of  ground 


Surface  impermeable  to  water  and  beat 


Figure  1.  Diagram  of  Columns  Arranged  in  a Row  Parallel  to 
the  Filtration  Flow. 


In  solving  the  problem,  we  mark  out  the  path  which  we  will  use: 
we  will  first  formulate  in  general  form  the  heat  balance  equation 
of  the  frozen  ground  wall,  we  then  determine  the  expressions  of 
its  components  and  substitute  them  into  the  initial  equation;  by 
integrating  the  differential  equation  which  is  obtained  in  this 
manner,  we  find  the  desired  solution  to  the  problem. 

3.  Heat  Balance  Equation  and  Its  Components 

a)  Heat  balance  equation.  The  intensity  of  the  heat  influx  from 
the  unfrozen  ground  varies  along  the  wall.  Therefore  the  thickness 
of  the  wall  along  its  length  is  different  (Figure  2).  Allowing  for 
the  fact  that  the  thickness  of  the  wall  changes  smoothly,  it  is 
possible  to  mentally  divide  the  wall  into  sections  of  finite 
length  and  to  assume  that  within  the  limit  of  each  section  it  has 
a uniform  thickness.  If  we  select  the  number  of  sections  to  be 
equal  to  the  number  of  columns,  then  the  length  of  each  section  will 
be  equal  to  S. 
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Figure  2.  Water-Permeab le  Ground  Frozen  by  a Row  of 
Columns  After  the  Cylinders  Solidify. 


The  heat  balance  equation  of  section  number  i of  the  frozen 
ground  has  the  following  form: 

Q..  = Q . - Q-.  . - Q (1) 

Hfx  xCX  xfl  X ^gx 

where  i is  the  number  of  the  section  (i  = 1,  2,  ...  I); 

Q^.  is  the  amount  of  heat  evolved  per  unit  of  time  in  freezing  water 
in  section  number  i; 

Qc^  is  the  amount  of  heat  which  is  extracted  per  unit  of  time 

from  the  ground  by  the  column  in  section  number  i (the  heat  absorbed 
by  the  freezihg  column) ; 

Qfi  ^ is  the  amount  of  heat  which  the  filtration  flow  gives  up  per 

unit  of  time  to  the  surface  of  the  frozen  ground  cylinder  in  section 
number  i; 

Q . is  the  amount  of  heat  which  is  extracted  per  unit  of  time 
gx  r 

from  the  frozen  ground  in  section  number  i (the  cooling  of  the  frozen 

ground  wall). 

Thermal  flow  Q . . as  will  be  shown  later,  can  be  taken  into  account 
gi 

by  simply  increasing  the  calculated  value  of  the  latent  ice  formation 
heat  which  belongs  to  the  expression  which  determines  the  thermal  flux 
. Because  of  this  fact  heat  balance  equation  (1)  acquires  the 

following  similar  form: 


Q..  = Q . - Q-.  . . (2 

vlx  ycx  yfx  X 

b)  Heat  given  off  at  the  freezing  boundary.  The  amount  of  heat 
which  is  given  off  per  unit  of  time  in  freezing  the  water  in  section  i 
is  equal  to 


where 


o is  the  amount  of  heat  given  off  in  freezing  a unit  of  ground 
volume; 

is  the  volume  of  frozen  ground  in  section  i; 
t is  time. 

Since 

Vt  — '2XjiS—-r~!t. 

1 1 0 

where  2X.  is  the  thickness  of  the  frozen  ground  wall  in  section  i,  we 
have 
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dV,  --2 hSdX,. 


(4) 


By  substituting  (4)  into  (3),  we  finally  find 

Qfl  — 2/tSz  . (5) 

d Z 

c)  The  heat  influx  from  the  unfrozen  ground.  The  heat  influx  from 
the  unfrozen  ground  in  section  i is  equal  to: 


where 

y^  is  the  distance  from  the  beginning  of  the  wall  to  the  beginning 
of  the  area  in  question; 

q~.  is  the  intensity  of  the  thermal  influx  from  the  unfrozen 
ground  to  the  wall  at  a distance  of  y from  the  beginning  of  the  wall. 

The  intensity  of  the  thermal  influx  can  be  found  in  the  following 
fashion. 

We  will  write  the  equation  of  the  established  heat  balance  of 
the  elementary  volume  of  unfrozen  ground,  and  we  will  then  convert 
this  equation  into  the  following  form: 


— = I i 6-1  , 

•>!<  c \ <>x-  ' • 

WV 


t is  the  temperature  of  the  ground; 

is  the  thermal  conductivity  factor  of  the  unfrozen  ground; 
c^  is  the  thermal  capacity  of  the  water; 

Yw  is  the  specific  weight  of  the  water; 
v is  the  speed  of  the  filtration  flow. 

Assuming  that  the  magnitude  of  the  thermal  flow  which  arises  due 
to  physical  heat  conductivity  and  is  directed  parallel  to  the  wall 
(along  axis  y)  is  considerably  less  than  the  value  of  the  heat  flow 
which  is  transferred  in  the  same  direction  by  the  water,  it  is 
possible  to  assume: 

££■ 


then  instead  of  (7)  we  will  have 


(3) 


r' 


where 


ww 


We  will  ignore  the  curvature  of  the  frozen  ground  wall  surface 
and  we  will  transfer  the  origin  of  the  coordinate  system  to  it; 
then 


when  x = 0 we  have  t = t^, 

when  x = ® we  have  t = T,  (9) 

when  y = 0 we  have  t = T. 

Equation  (8)  may  be  regarded  as  a one-dimensional  thermal 

conductivity  equation  in  which  the  time  variable  t is  replaced  by 
the  variable  coordinate  y.  The  solution  to  this  equation  under 
boundary  conditions  (9)  is  known  in  heat  transmission  theory  and 
looks  as  follows  [3,  p.  71;  4,  p.  91]: 


t — ( T—t0) 


(10) 


where  G is  the  integral  of  error  probability. 

From  the  last  expression  it  follows  that  the  temperature  gradient 
in  the  direction  of  axis  x is  equal  to 

.V*' 

ot  ___  r ~ - <„  'V 


At  the  surface  of  the  frozen  ground  wall  the  following  is  true 


<\i 

tlx 


j 
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Consequently, 
ground  to  the  wall 


the  intensity  of  the  thermal  influx  from  the  unfrozen 
will  be 


% 


'jAlsz'zl 

1 'V 


(11) 


Equation  (11)  was  formulated  in  1945  by  B.  V.  Proskuryakov  and 
was  used  by  him  in  solving  the  problem  of  freezing  water-permeable 
ground  [5], 


By  substituting  (11)  into  (6)  and  by  performing  integration,  we 
finally  find  the  desired  equation  for  determining  the  magnitude  of 
the  thermal  influx  from  the  unfrozen  ground  to  the  wall  in  section 
number  i: 


_ -J:f  i ( " — r'i) 

1‘  1 *6fi 


| • v • i-  S 


(12) 
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From  equation  (12)  it  follows  that  the  thermal  influx  from  the 
unfrozen  ground  does  not  change  during  freezing,  i.e.,  the  heat 
influx  is  a constant  value  for  every  section. 

d)  The  absorption  of  heat  by  freezing  columns.  The  heat  absorbed 
by  each  of  the  freezing  columns  in  the  row  depends  on  the  temperature 
of  the  brine  circulating  in  the  column  and  on  the  thermal  resistance 
of  frozen  ground  in  the  area  in  question.  Since  the  speed  with  which 
the  freezing  boundary  moves  is  extremely  slow,  the  temperature 
distribution  in  the  frozen  ground  is  quasistationary.  Therefore, 
in  the  general  case  we  have 


<?c, 


(13) 


where  is  the  thermal  resistance  of  the  frozen  ground  between  the 
column  and  the  freezing  boundary  in  section  number  i. 


By  utilizing  the  known  solution  to  the  hydromechanical  problem 
of  fluid  flow  to  an  infinite  row  of  ideal  holes  arranged  in  a 
band-like  stratum  [6,  p.  35],  we  can  write 


T,  - 


(14) 


where  n is  the  half-thickness  of  an  imaginary  wall  which  corresponds  to 
the  additional  thermal  resistance  which  arises  due  to  the  fact  that 
temperature  9 is  not  maintained  over  the  entire  surface  in  which  the 
row  of  columns  is  located,  but  only  on  the  surface  of  the  columns 
which  are  some  distance  apart: 


n - - 


(15) 


Aj  is  the  thermal  conductivity  factor  of  the  frozen  ground. 

Value  n is  fairly  considerable.  Thus  at  column  radius  rQ  = 0.05  m 

and  when  the  distance  between  the  columns  is  S = 1 m,  according  to 
(15)  we  find  that  n = 0.27  m.  Therefore  we  cannot  ignore  the  value  n. 

In  order  for  equation  (13)  which  expresses  the  heat  absorption 
of  the  column  to  be  single-valued  and  absolute,  it  is  necessary  to 
again  examine  the  value  of  brine  temperature  9. 

As  is  well  known,  the  temperature  of  the  brine  changes  during 
freezing.  These  changes  are  completely  regular  and  are  a function 
of  the  characteristics  of  the  refrigeration  device  and  the  increase 
during  freezing  in  the  thermal  resistance  of  the  frozen  ground.  Without 
making  a detailed  examination  here  of  the  problem  of  calculating  the 
course  of  the  brine  temperature,  we  will  simply  point  out  certain 
general  facts  in  this  area. 


The  operative  refrigeration  capacity  of  compression  refrigeration 
devices  primarily  depends  on  the  evaporation  temperature  of  the  coolant 


I 
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and  consequently  on  the  temperature  of  the  brine.  On  the  other  hand, 
the  thermal  absorption  of  the  columns  also  depends  on  the  brine 
temperature.  The  operative  refrigeration  capacity  of  the  devices 
should  always  be  equal  to  the  thermal  absorption  of  the  column. 

This  equality  occurs  (or  more  precisely  speaking,  arises)  at  a certain 
brine  temperature.  As  the  volume  of  the  frozen  ground  mass  increases 
and  consequent ly  as  the  thermal  resistance  of  the  frozen  ground 
increases,  this  equilibrium  is  disrupted  since  the  increased  thermal 
resistance  causes  a reduction  in  the  columns'  heat  absorption.  The 
production  of  cold  begins  to  exceed  the  requirement;  the  difference 
decreases  the  temperature  of  the  brine.  Equilibrium  is  restored  at 
a different  lower  brine  temperature.  This  is  how  we  find  the  dependence 
which  is  correct  when  the  columns  are  arranged  in  any  form: 

'■  -/VO.  (16) 

When  the  columns  are  arranged  in  a row,  after  the  frozen  ground 
cylinders  join  dependency  (16)  acquires  the  following  form: 

')  f(X\.  (17) 

The  reduction  in  brine  temperature  has  a limit  which  is  determined 
by  the  technical  operating  conditions  of  the  refrigeration  devices. 

When  this  limit,  which  corresponds  in  each  case  to  a certain  frozen 
ground  thermal  resistance  value,  is  reached,  subsequent  freezing  will 
take  place  at  a constant  minimum  brine  temperature: 


From  what  we  have  stated  it  follows  that  the  freezing  process  is 
broken  down  into  two  periods.  During  the  first,  initial  period  the 
temperature  of  the  brine  decreases,  and  during  the  second  period  the 
brine  temperature  is  constant. 

The  transition  from  the  first  period  to  the  second  is  determined 
by  substituting  (18)  into  (17): 


where  2Xf  is  the  value  of  the  wall  thickness  at  which  the  first 

period  of  freezing  stops  and  the  second  period  of  freezing  begins. 

It  must  be  noted  that  this  transition  to  the  second  period  of 
freezing  occurs  in  the  majority  of  cases  even  before  the  cylinders 
join,  and  therefore  freezing  after  joining  usually  occurs 
at  a constant  brine  temperature. 

e)  Consideration  of  the  quantity  of  heat  extracted  from  the  frozen 
ground  wall  when  it  is  cooled.  We  will  determine  the  quantity  of 
heat  which  is  extracted  per  unit  of  time  from  the  frozen  ground  in 
section  number  i when  the  wall  is  cooled.  In  doing  so  we  will  assume 
that  temperatures  9 is  not  maintained  on  the  column  surface,  but  rather 
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on  the  plane  in  which  the  row  of  columns  is  arranged.  In  examining 
expressions  (13)  and  (14)  which  determine  the  amount  of  heat  absorbed 
by  the  columns,  it  was  pointed  out  that  this  assumption  should  not 
be  adopted.  The  basis  for  introducing  this  assumption  in  this  case 
which  considerably  simplifies  the  mathematical  computations,  is 
the  following  considerations:  1)  the  value  of  component  Qoi  is 

usually  small  in  comparison  with  the  other  components  of  heat  balance 
equation  (1);  2)  this  assumption  obviously  leads  only  to  a small 
increase  in  the  calculated  value  of  Q ^ and  consequently  to  a very 

small  increase  in  the  calculated  value  of  freezing  time;  the  latter 
fact  provides  a certain  margin  in  designing  (the  time  calculated 
for  the  formation  of  a frozen  ground  wall  increases). 


Thermal  flux  Q ^ is  equal  to: 


g1 


dll', 


where  to\  is  the  heat  content  of  section  i of  the  frozen  ground  wall. 


Since 


IK  ->/ic(V(S  ,Y„ 


then  it  is  obvious  that 


Q . 

81 


(t  „-0) 


dX, 

a-  ' 


(19)  • 


where  c^  and  y^  are  the  heat  capacity  and  specific  weight  of  the  frozen 
ground . 

Now  let  us  compare  expression  (19)  with  equation  (5)  which  determines 
the  magnitude  of  thermal  flux  from  the  latent  ice  formation  heat.  It 
is  easy  to  see  that  if  we  stipulate 


then 


Qi 


o y/.sv 

'gl  -'O' 


(Ja- 


(21) 


This  is  how  we  can  allow  for  the  effect  of  thermal  flow  Q . by 

gi 

increasing  the  calculated  value  of  the  latent  ice  formation  heat  a. 

k.  Solving  the  Heat  8alance  Equation 

For  the  first  period  of  freezing  where  the  brine  temperature 
decreases,  heat  balance  equation  (1)  acquires  the  following  form 
when  che  expressions  of  its  components  (5)  and  (13)  are  substituted 
in  it,  with  allowance  for  (17)  and  (21): 
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//•V, 

dz 


(22) 


The  solution  to  the  stated  problem,  determining  the  time  required 

for  the  formation  of  a frozen  ground  wall,  is  found  by  integrating  the 

last  equation  within  the  limits  of  X.  = X.,  to  X.  = X._  and  x = x.  to 
^ 1 ll  1 i2  1 


T 


dX: 


(25) 


In  the  general  case  where  function  f(X^)  which  expresses  the 

brine  temperature  is  complex  in  form,  it  is  simplest  to  calculate 
the  value  of  the  integral  by  the  graphic  or  tabular  method. 


The  maximum  thickness  of  the  frozen  ground  wall  is  determined 
by  the  equation  which  expresses  the  equality  between  heat  absorption 
by  the  column  and  the  thermal  influx  from  unfrozen  ground: 


f n — / v X t ► 


- <V 


(24) 


In  cases  where  it  is  possible  to  assume  within  certain  limits  that 
the  column's  heat  absorption  is  constant,  i.e.,  that 

t o / f »X . i ^ 

— = Qa—  const. 

it  is  easy  to  find  from  (25)  that: 


Xr! 


^cii"  '-'fi.1 


-.)  Xtl. 


(25) 


We  will  now  derive  a computation  equation  for  the  second  freezing 
period  where  the  brine  temperature  is  constant.  Heat  balance  equation 
(1)  acquires  the  following  form: 


o/;cy ■ ■ ~'-1  — . ((°  ~ f'min)  n 
dz  A . - •:  fi‘ 


(2b) 


or,  by  separating  the  variables  we  find: 


JX„ 

rt'A  - :</> 


where 


."lAY 


(27) 


.hsi  i ll.,  — ‘i nun ' _ 

Hi-  ’ 


(23) 


By  integrating  the  latter  differential  equation  within  the  limits 

of  X.  = X.,  to  X.  = X.,  and  r = r,  to  r = r_,  we  obtain  the  following 

i ll  l :2  12 

calculation  equation: 
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(.v  — *)i: 


V - Xy 
' - .V  - 


-<*„  -V ,) 


(29) 


From  equation  (29)  it  is  evident  that  the  maximum  value  of  the 
half-thickness  of  the  wall  is  equal  to: 

■V,,  ; , .V.  (30) 

If  the  initial  value  is  > N,  then  obviously  formula  (29)  will 

indicate  the  melting  time  of  the  wall,  and  the  maximum  value  of  the 
half-width  of  the  wall,  as  before,  will  be  equal  to  N. 

In  order  to  make  computations  easier,  in  Figures  3 and  4 we  show 
computation  graphs  plotted  according  to  formula  (29) . 


Figure  3.  Graph  for  Determining  the  Time  Required  for  the 
Formation  of  a Frozen  Ground  Wall  (n  = 0.37  m) - 

5.  Replacing  a Row  of  Closed  Frozen  Ground  Cylinders  with  an  Equivalent 
Wal  I 

In  the  computation  equations  derived  in  the  previous  paragraph, 
the  initial  wall  thickness  was  arbitrary.  We  are  interested  in 
the  speed  of  freezing  after  the  time  when  the  frozen  ground  cylinders 
joined.  Therefore  we  need  to  determine  the  wall  thickness  during 
this  time.  We  will  assume  that  instead  of  the  actual  wall  which  is 
composed  of  cylinders,  we  are  dealing  with  a flat  rectangular  cross- 
-section  wall.  We  will  determine  what  thickness  of  this  "equivalent" 
wall  should  be.  In  doing  so  it  is  necessary  to  observe  two  conditions 
which  consist  of  the  fact  that  the  volume  of  frozen  ground  and  the 
thermal  resistance  between  the  column  and  the  surface  of  the  frozen 
ground  must  remain  constant. 
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(32) 


2X»  — - ~ S - 0.70S. 


We  will  write  an  equation  which  corresponds  to  the 
The  thermal  resistance  of  the  frozen  ground  cylinder  ^ 


second  condition, 
is  equal  to 


consequently,  at  the  moment  the  cylinders  join 


(33) 


The  thermal  resistance  of  the  equivalent  wall  with  a length 

which  is  equal  to  the  distance  between  the  column  axes  S can  be  presented 
in  the  following  form  in  accordance  with  equation  (14): 


T. 


(34) 


In  accordance  with  the  second  condition  1*  = f . 

c r 

By  substituting  expressions  (33)  and  (34)  into  this  equation, 
we  find 


-r0 


or,  keeping  in  mind  the  fact  that 


we  obtain 


from  this  we  finally  find  that  in  accordance  with  the  second  condition 
the  thickness  of  the  equivalent  wall  may  be  equal  to: 


A comparison  of  expressions  (32)  and  (35)  shows  that  they  provide 
results  which  diverge  by  less  than  8%.  Therefore,  by  satisfying  the 
above-mentioned  two  conditions  fairly  precisely,  we  can  assume: 


A'„  =•-  0.  ;SS.  (36> 

From  this  expression  it  follows  that  the  thickness  of  an  "equivalent" 
wall  is  approximately  equal  to  three-fourths  of  the  distance  between  the 
axes  of  the  freezing  columns. 
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6.  Sample  Calculation 


We  will  determine  the  time  it  takes  for  a frozen  ground  wall 
to  form  beginning  at  the  time  when  the  cylinders  join  and  before 
the  design  thickness,  equal  to  two  meters,  is  reached. 

The  following  are  known: 


speed  of  filtration  flow 

temperature  of  filtration  flow., 
freezing  temperature  of  ground.. 

specific  weight  of  frozen  ground 

thermal  conductivity  factor  of 
frozen  ground 

heat  capacity  of  frozen  ground.. 

thermal  conductivity  factor  of 
unfrozen  ground 

latent  ice  formation  per  unit  of 

ground  volume 

heat  capacity  of  water 

specific  weight  of  water 

radius  of  freezing  column 

length  of  one  freezing  column... 
distance  between  column  axes.... 
number  of  freezing  columns 


v = 0.0S  m/hr 
T = 4°C 

C0  = °°C 

Aj  = 1600  kg/mJ 

Aj  = kcal/m-hr’deg 
Cj  = 0.35  kcal/kg-deg 

A?  = 1 kcal/m-hr'deg 

0 = 23,000  kcal/m-3 
c^  = 1 kcal/kg-deg 

A =1  kg/m’’ 

rQ  = 0.05  m 

h = 10  m 
S = 1 m 

1 = 25. 


Freezing  occurs  during  the  second  period  when  the  brine  temperature 

is  constant  0 . = -20°C. 

min 

Let  us  primarily  examine  the  first  section  (i  = 1). 


According  to  (12)  the  heat  influx  from  the  unfrozen  ground  is 
equal  to: 


Ofn- 


IftAj  [ T — 


■a- 


ll(Mil  I — h!()  Jkpal 
i hr 

r "i.K'Oh.on' 


According  to  expression  (15)  we  find: 

: :n  - o.  ,-  , 

• r 

Now  from  (28)  we  can  determine  the  value  of  N: 

,V  = zLVHl-L?’.  (M7  , , 

V til  r,'!0 

The  maximum  wall  thickness  is  equal  to  2X  j. . =.  2.V  --  1.7'.  , i.e., 

less  than  expected. 
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rT*^fr' 


•»»  • 


Let  us  examine  the  fifth  section  (i  = 5).  In  the  same  way  as 
we  calculated  the  first  section,  we  find 

4-ltM-4  (I/4TT  - 1 ~,)=  io4 
- 1 ’ 

" 1 1000  0,05 


151 


Here  in  the  fifth  section  the  maximum  thickness 

already  considerably  exceeds  the  expected  value  2X  = 2 m.  From  formula 
(29)  we  will  determine  the  necessary  time  of  active  freezing.  However, 
we  first  note  that  according  to  expression  (36)  the  initial  wall  thickness 
is  equal  to  2X^  = 0.76  S = 0.76  m;  the  value  of  a'  is  determined  from 

equation  (20) 

s'  = - £■,•(,  — — = 23000  - 0.35-1600  — = 26600  kcai/ra3. 

o 1 


According  to  expression  (27)  the  value  of  B is  equal  to: 

B - = Lil 2 1,9. lo~4  m/hr 

2hS-.’  2-10- 1-2660*)  - 

Keeping  in  mind  the  fact  that  Xj  = Xq  = 0.55  m and  X2  = 1 m,  according 
to  formula  (29)  we  find: 

= L . f ( |,S3  _ 0.37)  In  — - (1  - 0.38)1  ~ 595  hrs  - 25  days  . 

2,69-10- '[  4.83-1  j 


Let  us  examine  the  last  section  of  the  wall  (i 
= 040  ( \ 24  -T-T  - I/Jt)  = 64  kcal/hr  ; 


II 

Is 

1 

© 

CO 

= 12,13 

64 

1,4  = 1.12 

• 10  ~4  m/hi; 

372000 

25).  We  find: 


IQiJOO 

1,12 


[(12,13  — 0.37)  In 
I.  12.13-1 


- (1  - 0,38) 


495  hrs  r 21  days  . 


This  calculation  shows  that  the  difficulty  consists  of  forming 
a frozen  ground  wall  in  the  first  (initial)  section.  It  is  obvious 
that  here  special  measures  should  be  undertaken,  for  instance,  the 
implantation  of  several  additional  freezing  columns.  As  a whole, 
however,  the  formation  of  a frozen  ground  wall  after  joining 
of  the  cylinders  will  take  approximately  25  days. 
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7.  Conclusions 


The  main  results  of  this  article  amount  to  the  following. 

1.  Calculated  dependencies  are  obtained  which  make  it  possible 
to  determine  the  freezing  and  thawing  time  of  a frozen  ground  wall 
which,  after  frozen  ground  cylinders  are  formed,  is  located  parallel 
to  the  motion  of  the  filtration  flow: 

a)  for  the  first  freezing  period  when  the  brine  temperature 
changes,  computation  equation  (23);  in  the  special  case  where  the 
heat  absorption  of  the  columns  is  constant,  computation  equation  (25); 

b)  for  the  second  freezing  period  where  the  brine  temperature 
is  constant,  computation  equation  (29). 

These  computation  equations  differ  from  similar  formulas  mainly 
due  to  the  fact  that  they  make  greater  allowance  for  the  values  of 
thermal  influx  from  the  unfrozen  ground  and  the  heat  absorbed  by  the 
columns. 

2.  Computation  equations  (24)  and  (30)  determine  the  maximum 

value  of  the  wall  thickness  towards  which  it  will  tend  when  freezing 

(or  thawing)  takes  an  infinite  amount  of  time. 

3.  The  intensity  of  the  thermal  influx  from  unfrozen  water-permeable 
ground  varies  along  the  length  of  the  wall;  as  water  moves  along  the 
wall,  the  water  is  cooled  and  the  intensity  of  the  thermal  influx 
decreases.  From  this  it  follows,  in  the  first  place,  that  it  is 

most  difficult  to  create  the  initial  (with  regard  to  the  flow)  section 
of  the  frozen  ground  wall,  and  in  the  second  place  that  the  creation 
of  the  wall  must  begin  with  the  upper  sections  since  this  makes  it 
easier  to  create  the  lower  sections. 

4.  The  value  of  quantity  .4  - )...c  — f.,)  may  make  it 

possible  to  make  a preliminary  evaluation  of  the  difficulty  of  freezing 

water-permeable  soils.  The  greater  the  value  of  A,  the  more  difficult 
it  is  for  freezing  to  take  place.  At  values  A = 100-120 

kcal/m^” • hr  and  above,  freezing  is  associated  with  considerable 
difficulties  and  therefore  may  prove  to  be  economically  unfeasible, 
as  well  as  sometimes  even  technically  impossible. 

In  conclusion  I would  like  to  note  that  in  writing  this  article 
Professor  (Doctor  of  Technical  Science)  B.  V.  Proskurvakov  provided 
many  valuable  comments. 
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